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Fermions play a special role in homogeneous models of quantum cosmology because the exclusion 
principle prevents them from forming sizable matter contributions. They can thus describe the mat- 
ter ingredients only truly microscopically and it is not possible to avoid strong quantum regimes by 
positing a large matter content. Moreover, possible parity violating effects are important especially 
in loop quantum cosmology whose basic object is a difference equation for the wave function of the 
universe defined on a discrete space of triads. The two orientations of a triad are interchanged by a 
parity transformation, which leaves the difference equation invariant for ordinary matter. Here, we 
revisit and extend loop quantum cosmology by introducing fermions and the gravitational torsion 
CNl ' they imply, which renders the parity issue non-trivial. A treatable locally rotationally symmetric 

Bianchi model is introduced which clearly shows the role of parity. General wave functions cannot 
be parity-even or odd, and parity violating effects in matter influence the microscopic big bang 
■ transition which replaces the classical singularity in loop quantum cosmology. 

l> : 

I. INTRODUCTION 

o ! 

q-( Most cosmological models — classical or quantum — introduce the matter ingredients of the universe as bosonic 
fields, in particular scalar ones. While this provides a good measure for the implications of matter energy on space- 
^q' time, some effects of realistic fermionic particles may be overlooked. Especially in homogeneous models of quantum 
i— i, cosmology there is an important difference between bosonic and fermionic models: the exclusion principle forbids large 
matter energies when symmetry reduction leaves only a few, finitely many fermionic degrees of freedom. A massive 
universe can then be obtained by only two possibilities: a homogeneous description with many different fermionic 
I ■ species, or inhomogeneity with many local degrees of freedom of a few species (as in [![). 
^v^j I Both options differ from what is modeled by large values of homogeneous bosonic fields which rather resemble a 
Bose-Einstein condensate of many identical excitations. In fact, fermion condensates have been suggested for such a 
purpose, with characteristic effects 2. 3]. This is an effective picture starting from an inhomogeneous perspective in 
which fermions condense under certain conditions, after which a symmetry reduction can be done. It differs from a 
fundamental description from fermions in quantum cosmology where constraints due to the exclusion principle cannot 
as easily be avoided. Potentially fundamental mechanisms which rely on a large amount of bosonic matter, such as 
bounce scenarios to avoid the big bang singularity, have to be reanalyzed if matter is fermionic. A truly microscopic 
description will then be achieved. Here, we perform an analysis of the role of fermions in loop quantum cosmology. 

Loop quantum cosmology [H provides a general mechanism for fundamental singularity resolution Q . Commonly 
in quantum cosmology, when volume is used as an intrinsic measure of time, evolution must stop at the classical 
singularity where the volume vanishes. In loop quantum cosmology, by contrast, the timeline is naturally extended, 
first at the kinematical quantum level, by including orientation into the basic variables: The (densitized) triad knows 
about the size as well as the orientation of the universe which make it take all real values, not just positive ones. 
Vanishing size is then no longer a boundary but an interior point of minisuperspace. 

What is more, even dynamically the classical singularity is removed because the quantized Hamiltonian constraint 
equation uniquely extends any wave function defined on minisuperspace across the subset of vanishing sizes. Dynamics 
is dictated by a difference equation which remains regular where classical relativity and Wheeler-DeWitt quantum 
cosmology would reach their limits. Rather than being singular, the big bang transition then appears as a place where 
space flips its orientation — turning its inside out — while it changes from being contracting to being expanding. 1 
Still, the region of vanishing volume does remain special in the underlying recurrence scheme. Some coefficients of the 



Sometimes it is suggested that this small-volume regime is avoided altogether because wave packets may turn around in a bounce at 
some minimal non-zero volume. This indeed happens for homogeneous models containing sufficiently much kinetic energy of matter 
1 1 Of) - However, this does not appear as a general mechanism which would be valid in this form for generic quantum states or for 
inhomogeneous situations. Fundamental singularity resolution which deals with the wave function right at vanishing volume is thus 
required. 
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difference equation can vanish at labels corresponding to zero volume, which leads to consistency conditions implied 
by the dynamical law [ill) E3] • This is welcome because, at least partially, it frees one from having to pose initial 
values for a wave function independently of the dynamics. The wave function of the universe is restricted by the 
theory alone, relaxing the need to pick one solution among many which could correspond to our universe. 

While the set of configurations of vanishing volume is not a boundary within the theory, in the presence of ordinary 
matter one may choose to consider parity transformations as large gauge transformations which complete the gauge 
group of triad rotations to all orthogonal transformations. Then, one would restrict solutions to only those states 
which are either even or odd under parity reversal. This would essentially factor out the orientation degree of freedom 
introduced by the use of triad variables, and again demote the set of vanishing sizes to a boundary rather than an 
interior regime. This factoring has indeed been assumed in recent constructions of physical Hilbert spaces for specific 
isotropic models and the corresponding intuitive bounce pictures based on (l3| . 

But if this is used crucially for the constructions, what happens if more realistic matter is included which, as we 
know from particle physics, cannot be parity invariant? Do properties of the specific solutions based on the assumption 
of reflection symmetry depend on the conservation of parity by matter, and if so, how reliable are the conclusions 
drawn from this assumption? Only the inclusion of parity violating terms, at least as a possibility, can provide a 
sufficiently general mechanism of singularity resolution. 

It may also give rise to new effects related to the role of parity violation in the big bang transition. If this were 
to happen, an intriguing new link between particle physics and quantum gravity would result. Seeing whether this 
is indeed the case requires the introduction of fermions, which is available in loop quantum gravity [lj, HH [l(| (see 
also [lj], EM EH). In general, however, the parity behavior of loop quantum gravity is highly non-trivial due to the 
fact that the basic variable conjugate to the densitized triad, namely the Ashtekar-Barbero connection, is the sum 
of a parity-even and a parity-odd term. It does not have a simple parity behavior and, moreover, it appears in 
quantized expressions only non-linearly through holonomies. Even in vacuum, this makes a direct demonstration of 
parity invariance of loop quantum gravity — or the lack thereof — very complicated [la ] . 

In this article, we introduce a homogeneous model which allows one to analyze the parity behavior in a clear-cut way. 
At the same time, the model is amenable to the techniques which have been proven useful for explicit constructions 
of Hamiltonian constraint equations through the difference equations of loop quantum cosmology 0, [M The 
microscopic nature of fermions due to the exclusion principle is explicitly realized. As we will see, quantization of this 
model does not introduce unexpected parity violations in the absence of classical parity violations. But the inclusion 
of parity violating matter interactions is possible, which can be used to illustrate the role of parity for singularity 
removal. Then indeed, wave functions change under triad reflections. The big bang transition through vanishing sizes 
is a non-trivial event, which represents true local evolution in internal time rather than merely the application of a 
symmetry transformation. 

II. CLASSICAL SYMMETRY REDUCTION 

In this section, we provide the formulation of symmetry reduced cosmological models which may have torsion due 
to the presence of fermions. We follow the symmetry reduction of torsion- free Bianchi class A models 0, [U , 
combined with the canonical formulation of gravity with fermions [l4l fl5L [l6l . [23| ; the general formulation without 
symmetry is summarized in the Appendix. Here, we combine these research lines and explore the symmetry reduction 
of gravity coupled to fermions in a first-order formalism, implying a theory with torsion. As we will see, there are 
non-trivial changes in the underlying equations, such that the analysis done here provides a crucial consistency test 
of the robustness of existing models. At the same time, we clarify the constructions of loop quantum cosmology Q 
from the viewpoint of some recent developments. 

A. Diagonalization 

Bianchi class A models constitute all homogeneous models with a symmetry group S acting freely on the space 
manifold £ = S and for which standard Hamiltonian formulations exist. The symmetry group is characterized by its 
structure constants Cj K , which for class A models satisfy Cjj — [24[ and can be parameterized as Cfj — e K I J rS K ' 1 
with three coefficients n 1 which either vanish or take values ±1. Some of these models can be reduced further by 
imposing rotational symmetry with one axis (where S has isotropy group F = U(l)) or even isotropy (F = SO(3)). 
Later in this paper we will present a locally rotationally symmetric (LRS |25| ) model with torsion in detail. 

The action of a symmetry group 5onS provides invariant 1-forms u>^ which are used for the reduction of Ashtekar- 
Barbero variables. For each s € S, they satisfy s*uj^ = Ad(s)jOj^ or, in terms of the Lie-algebra valued 1-form 
il a := lu^Tj with generators Tj of S, s*f2 a = s _1 f2 a s. The left invariant 1-forms then yield the decomposition 
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A l a = V 1 ^4>\io I a of an invariant connection with spatially constant coefficients ifi] (see the Appendix of Q for more 
details on invariant connections). Here, we have explicitly included a factor of Vq = J d 3 x\ det(w^)| of the spatial 
coordinate volume (or the volume of any finite region used to define the homogeneous variables) as it will be convenient 

later on. A corresponding decomposition of the densitized triad is given by Ef = V~ 2/3 p{Xf with Xf being densitized 
left invariant vector fields dual to the 1-forms: ui^Xj = 5j \ det(u!^)\. The symplectic structure of the reduced model 
is given by 

7*3# (1) 

as it follows from (7/4) _1 J d 3 xA i a Ef — (•yK)~ 1 (f) i I pf. 

For the purpose of loop quantization, it is useful to further reduce the number of independent components of 
the invariant connection and its conjugate momentum. In some cases, this will allow very explicit calculations of 
matrix elements of the Hamiltonian constraint and the difference equation it implies for physical states [20I ] . Both the 
connection and the densitized triad can be cast into diagonal form 

4 = V V V)A^f , E? = V 2 /y*)Af X% (2) 

with six spatially constant coefficients ci and p 1 which are considered as the only dynamical components while 
A € SO(3) is fixed up to gauge transformations. Using the same A} for A\ and Ef is consistent with the Gauss 
constraint for diagonal torsion-free Bianchi class A models which is then solved identically From the diagonal 
densitized triad, moreover, we find the co-triad e l a = V Q l ' 3 aiK) A 1 k lo^ with \a% \ = \/\p 2 p 3 /p 1 \ and cyclic. It determines 
the diagonal anisotropic spatial metric 

q a b = e l a e\ = V^^afrfucjlujg = qu^W h 

with three independent scale factors V Q \aj\. 

By construction, c/, p 1 and aj are independent of coordinates as long as the diagonalized homogeneous form is 
respected. In particular in a Bianchi I model where lu^ — d a x T = 5^ in terms of Cartesian coordinates x 1 , spatial 
coordinates can be rescaled arbitrarily without affecting the basic variables. However, the specific values do depend 
on Vq and the choice of the integration volume. Obviously, the Vb-dependence is a consequence of the symmetry 
reduction to homogeneity, since Vq does not occur at all in an inhomogeneous framework. Thus, the dependence 
has to be interpreted with care especially after quantization where, fundamentally, the relation to coordinates is lost. 
As a consequence, the role of Vq cannot be properly understood if considerations are limited to purely homogeneous 
models because only the reduction from inhomogeneity shows how Vq enters; see [26[ for a discussion from the point 
of view of inhomogeneous states. 

Note that p 1 and ax are allowed to take negative values to represent different triad orientations while the orientation 
of A e SO (3) is fixed. A parity transformation then simply implies p 1 1— > — p 1 for the triad components (leaving 
coordinates unchanged), while the transformation of the c/ is in general more complicated. In fact, we have A l a = 
T l a + jKq with the parity-even torsion-free spin connection 

f^i6^e fc fe (29 [b < ]+ej c el9 h ei) (3) 

and the odd extrinsic curvature K % a = K a f,e\. In the torsion-free case, it follows from ([3]) that the homogeneous spin 
connection can be expressed as T l a — ^ik) 

A^wf HI with 

Tj = - \^-n' J H — —n K — n 1 ] for indices such that £ijk = 1 (4) 

2 \a,K aj ajaK J 

and the same A} as used for the densitized triad. Similarly K % a = KmAj u> l a also with the same A}. Then, cj = Tj+jKi 
does not have a straightforward parity behavior unless Ti — (as in the Bianchi I model). 

The diagonalization is sufficient to capture the crucial dynamical behavior of Bianchi models, such as the approach to 
a singularity. For the quantization, it has the advantage that it reduces SU(2) to U(l) 3 : holonomies of a homogeneous 

connection, computed along curves generated by the invariant vector fields Xf, take the form — exp(/z</>}Ti) with 
a real number /i depending, e.g., on the coordinate length of a curve used to compute the holonomy. For (f> l j — C(/)Aj, 
we have 



hf'> = exp(^/C(/)AjTi) = cos (|At/C( 7) ) + 2A}r; sin {\^iC(d) 



(5) 
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While any SU(2)-holonomy along Xf can be written in this way, 2 the diagonalization implies that A\ becomes a 
mere background quantity not subject to dynamics. Thus, it is sufficient to consider only the simple commuting 
exponentials exp(z^/C(/)) to separate diagonal connections. After a loop quantization, as we will see in detail below, 
this will have the implication that a triad representation exists, which simplifies the analysis of dynamics considerably. 
In fact, triad operators will simply be p 1 = —i"/£pd/dci, with the Planck length £p = \/kK, which form a complete 
commuting set. Their eigenstates 

(ci,c 2 ,c 3 |/ii,M2,M3) = exp(ii(/iici + /i 2 c 2 + pt 3 c 3 )) 
(written here in the connection representation) form an orthonormal basis such that the coefficients in 

form the triad representation of arbitrary states. This explicit representation, which becomes available only after 
diagonalization |2(j, has been the basis of all investigations so far in homogeneous loop quantum cosmology. As we 
will see in this article, arriving at such a representation is less trivial in the presence of torsion. 

B. Torsion effects 

This scheme of diagonalization of the basic torsion-free gravitational variables relies on the fact that both the 
connection and its conjugate momentum can be diagonalizcd with the same Aj. In other words, the su(2) valued 
connection and its conjugate momentum are parallel to each other in the tangent space of the internal symmetry 
group. This can be seen from the torsion-free Gauss constraint which expressed in terms of the diagonalized variables 
takes the form p^c^eyfcAjA^ = and is identically satisfied. However, the presence of torsion via the axial fermion 
current J;, as summarized in the Appendix, enters the Gauss constraint (IA7p implying that 

4pUm = ^\det{p I j )\J i . (6) 

For <ft\ — C(/)A} and p\ = p^Aj as above, this would only allow vanishing spatial components of the fermion current 
and severely restrict the allowed models. This situation becomes more obvious if we try to express the spin connection 
including its torsion contribution as T l a = T^ K )A 1 k uj^ with the same Aj as used for the triad: One can easily verify 
that the partial torsion contribution (|A6j) to the connection cannot be expressed as C* = C^k)A 1 k oj^ if J 1 ^ 0. Then 
also the Ashtekar-Barbero connection cannot be diagonal in the same basis. Therefore, our first result is that the 
presence of torsion does not allow us to diagonalize both canonical variables, i.e. the connection and the densitized 
triad, simultaneously. 

Moreover, fermion terms require us to use a connection A l a in (I Al 1|) which carries an extra term compared to the 
Ashtekar-Barbero connection, depending on the fermion current. We then write the new diagonal variables as 

^ = ^ 17 V)A>f , E? = V- 2/ V K) TFXk (?) 

where in general T/ ^ A\. Not both A\ and T* can be fixed because partially they are determined by dynamical 
fields as, e.g., per the Gauss constraint ©. This has an immediate implication for the symplectic structure because 
cj and p 1 will no longer be canonically conjugate: 

/ d 3 ^£,4 = V ^TlL t ( C(0 Aj) = p«A (c (0 A*T/) - c {l) p^A)C t T( . (8) 
Js 

Thus, it is not cj which is conjugate to p 1 but cm A}T/ / ' ) . This is not a pure connection component but depends on the 
relative angles between the connection direction A^ and the triad direction T? in internal space. (It is not possible to 



2 General curves do not provide this simple form. For instance, along X™ + Xj holonomies are not of the (almost) periodic form in 
ci or C2 (but in y t- 2 ~ If curves are considered which are not even straight with respect to the given symmetry, the behavior is 
more complicated due to path ordering and do not give rise to almost periodic functions (StJ. However, such curves do not play a role 
in the kinematical symmetry reduction, which uses the given set of Xf to introduce particular quantum geometries, just like classical 
symmetric metrics which are used in adapted coordinates but can look complicated in arbitrary coordinates. 
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fix both A} and T/ because this would require six parameters while the Gauss constraint allows one to fix only three.) 
Moreover, some of the angles enter the symplectic structure as independent variables. We can, for instance, (Euler) 
parameterize T- as the matrix T((f>j) = exp((/> 3 T 3 ) exp(0 2 Ti) exp(</>iT 3 ) using generators T[ of SO(3). Inserting this 
in ([5]) shows that the angles (f>i acquire canonical momenta given in terms of the angles in Aj, e.g. 4>\ being conjugate 
to — tr((c • A)(p • T(tfii 4- 7r/2, cf> 2 , 4>3))), where c and p here denote the diagonal matrices with components cj and p 1 , 
respectively. (Taking a derivative of T{<f>j) amounts to switching sines and cosines, which is the same as shifting an 
angle by n/2.) 

The corresponding phase space and the constrained system defined on it is rather involved, and so we consider a 
more special case which still allows the non-trivial implications of torsion to be seen: We are interested in the case 
where the presence of a fermion current is the sole reason for anisotropy, while the 2-dimensional space transversal 
to the spatial current is rotationally invariant. We can then assume that there are bases for A\ and Ef, respectively, 
such that 

/ 1 \ / 1 \ 

A j j = cosp -ship , Tf — cos0 sin0 , (9) 
y ship cos/3 / y — sm(j) cos0 / 

where p and </> are the only non-vanishing rotation angles. As we will demonstrate below, this allows non-trivial 
solutions where the fermion current is aligned in the 1-direction. The Liouville term in the action can then be 
expressed as 



l - [d?xEfC t A\ = i-pWA(c (/) A^)- W «A}A^ 

— (tip 1 + C t (c 2 cos(p - 4>))p 2 + C t (c 3 cos(p - 0))p 3 - <fi(c 2 p 2 + c 3 p 3 )sin(p - <f>)) 

— (cip 1 + c 2 p 2 + d 3 p 3 + , (10) 



7« Jt, 7 k 
1 

7K 



where we introduced 

c 2 = c 2 cos(p - (f>) , c 3 = c 3 cos(p - (f>) , = -{c 2 p 2 + c 3 p 3 )sin(p - (f>) . (11) 

In these components, the symplectic structure is 

{ci,p 1 }=7K , {c 2 ,p 2 }=7K , {c 3 ,p 3 }=7K , {0,p }= 7 K. (12) 

Notice that the presence of torsion at this stage introduces a new kinematical degree of freedom (p. It will be removed 
after solving the Gauss constraint ([5]), which is now non-trivial. 

There is a useful interpretation of the canonical variables in the presence of torsion: We can write, e.g., 

~c 2 = c 2 cos(p -(/>) = c 2 AlT 2 = 4>\T 2 

in terms of the general homogeneous coefficients (f>\ = C(/)A}. Since T/ gives the direction of Ef , we can interpret c 2 
as a component 

V c 2 - A a E t 

V p 2 

of the projection of A\ onto Ef. In the absence of torsion, this would be a pure connection component because A\ 
and Ef would be parallel. With torsion, however, c 2 is only part of an ^-component: Using the expression fA13]) . the 
projection removes the term e l kie^J l perpendicular to Ef which happens to be the torsion contribution to extrinsic 
curvature. Moreover, the projection transversal to Ef is just (half of) the variable p^ due to the sine, which thus takes 
a value equal to the torsion contribution. This agrees with the solution of the Gauss constraint (fH)) below. Recall that 
the identification of the torsion contribution to extrinsic curvature used in (j A13|) cannot be completed without partially 
solving equations of motion. In the projection defining c 2 and c 3 , on the other hand, no equations of motion have been 
used. Thus, these canonical variables which we are naturally led to at the basic kinematical level present torsion- free 
contributions without explicitly splitting off torsion. (Something similar happens in inhomogeneous models such as 
spherical symmetry 28] or Gowdy models 29] . There it is spin connection contributions that are split off by a natural 
definition of canonical variables which then allows a manageable loop quantization.) 
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C. Reduced constraints 



In terms of the diagonal variables the Gauss constraint (|A7[) becomes 



Gi = —e iik c {I)V ^k\Ti - iybWj^J, = -^i(c 2 p 2 + c 3 p 3 )sin(p - <)>) - VlP ? P ' ifjj = . (13) 
7^£ a 7^ £ 

For z = 2, 3, it thus implies J 2 = = J 3 while the remaining condition 
relates J\ to p^: 

= y V\pW\Ji =■■ , (14) 

where Ji — + X a'iX denotes the densitized axial fermion current (which is bilinear in half-densitized fermions 

£ and x)- With the choice © of bases the fermion current Ji is aligned along the first (fixed) internal direction: 
J 2 — J 3 — 0. This defines a specific class of models with a non-trivial spatial fermion current, as J\ may be non-zero. 
Similarly, the diffeomorphism constraint (|A8|) can be written as 

V a N a = -cfj^pj N 1 = N 1 (n 2 c 2 p 2 + n 3 c 3 p 3 )sm(p - <j>) = , (15) 

where iV a = N 1 Xf with N 1 constant and Cfj = e K IJ n^ K ^ to specify different Bianchi class A models are used. 
We have also imposed that the partial derivatives of spinor fields vanish in a homogeneous model, e.g. d a tp = 0. 
A conclusion to be drawn from (TT3]) and (|15[) is that torsion is strongly restricted in Bianchi Class A models with 
n 2 + n 3 7^ since this implies that p^ — jkJ'i/2 — 0, and thus all spatial components of the axial vector current 
vanish. 

Finally, the Hamiltonian constraint (|A9[) is 



Bianchi 



yj\p 1 p 2 p 3 



(n 1 cip 2 p 3 + n 2 c 2 p 1 p 3 cos(p — 0) + n 3 c 3 p 2 p 1 cos(p — 0)) 



VbWi 

k-^- 2 (i + 1 2 ) 



C2P 2 C 3 P 3 ) 



\J\P X P 2 P 3 
1 



(cip 1 c 2 p 2 cos(p — 4>) + cip 1 c 3 p 3 cos{p - i 

(ci - f 1 )p\f 2 p 2 + f 3P 3 )cos(p - 0) - (ca - T 2 ) P 2 f 3 p 3 



^7(c 2 p 2 +c 3 p 3 ; 



sin(p - ^ + (t.p 1 + (t 2 p 2 + f 3 p 3 )cos(p - 0)) J° 



(n 1 


p 2 p 3 




p x p 3 




2 1 
p^p 1 






P 1 


+ n 2 


p2 


+ n 3 


p3 


) 



J° 



3jK0 



~0 



16v/bWI(l+7 2 ) 



2 7 £ 3 - - + 2 7 ' 
a 



16 v / \p 1 p 2 p 3 \ \a I + 7 
2 ) Jl 



J 2 



(16) 



where a, /? and are defined in the Appendix. It is important to emphasize that since Y % a is not diagonalized in either 
Aj or T/ in the presence of torsion, the Hamiltonian constraint in (| A9I) expressed in terms of F^ by splitting torsion 
from the spin connection is essential to obtain a controlled loop quantization as will be shown below. 



D. The Bianchi I LRS Model with Torsion 



If there is an isotropy group F = U(l) for the action of the symmetry group 5, one obtains locally rotationally 
symmetric (LRS) models. Therefore, two of the diagonal components of the connection as well as of the triad, e.g. 
the second two for definitiveness, have to equal each other and only two degrees of freedom are left which we choose 
to be (ci,p x ) and (c 2 ,p 2 ) embedded into the general Bianchi model by 



(ci,c 2 ) (ci,c 2 ,c 3 ) = (ci,c 2 ,c 2 ) , (p\p 2 ) b\p 2 ,P 3 ) = b\p 2 ,P 2 ) 
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The symplectic structure can be pulled back by this embedding providing Poisson brackets 

{ci,p 1 }=7K, {c 2 ,p 2 } = ^7K, {(f),p<f,} =JK 



(17) 



-2c, 



p 2 siri(p — if)). (Solutions of this symmetry type in the presence of torsion due 



from (TT2]) , where p$ is now p^ 
to spin fluids have been studied in [30|, |31|.) 

For the LRS model, the diffcomorphism and the Hamiltonian constraints, (fT"5|) and (fl"6)) respectively, further reduce 

to 



and 



V a N a = -^-N^n 2 + n 3 ) PrP = 



(18) 



LRS 



/ 1 ( 2\ z i 2-12, 3-21 -2 / o 1- 2 , (~ 1\ A , 1 2 

= (n ci(p ) +n c 2 p p + n°c 2 p p - 7 I 2cip c 2 p + (c 2 p ) + 



K -V 2 (l+7 2 ) 



,2 1 



2|£ 2 P 21 



(d-ro/^ + ra^) 



sJpl + ±{c2P 2 ) 2 



2|pa|vVI 



7P^i - A fV + (f 2 p 2 + f 3 ;p 3 



2|£ 2 P 2 



^/p 2 +4(5 2 p 2 ) 2 



'sgn(c 2 p 2 ) 



^+4(c 2 p 2 ) 2 -r 2 ^ \T 3 p 



7 



; 2 7 /3 ( 3 - -i- + 2 7 2 

where we have used the definitions of c 2 and p* to write 



2 ^ -r2 



(19) 



cos(p — <j)) 



2|5 2P 21 



^P% + 4(5 2 p 2 



(20) 



To allow a non-vanishing J7i and to be specific, we work from now on with the Bianchi I model. Here, the 
diffcomorphism constraint (p~8|) vanishes identically and does not impose any restriction on p$. This has the additional 



advantage that the resulting Hamiltonian constraint will be free of terms such as ypfj, + 4(c 2 p 2 ) 2 , which lack simple 

quantizations. (While there are well-defined operators with this classical limit, given that both p^ and (c 2 p 2 ) 2 would 
be mutually commuting positive operators whose square root can be taken after summing them, not all the operators 
involved have discrete spectra. Thus, it would not be straightforward to compute explicit matrix elements of the 
square root operator which would be required for the quantized Hamiltonian. Once the square root is quantized, its 

inverse in (|2H)) could easily be obtained from 2~/np 2 cos(p ~ 4>) — {^/p^ + 4(c 2 p 2 ) 2 ,p 2 }.) 

For the Bianchi I LRS model, we then have V 1 = and thus the Hamiltonian constraint is finally given by 



I LRS 



\p 2 \VW\ 

37K0 



2c 1 p 1 c 2 p 2 + (c 2 p 



7, „2\ 



7 



2b 2 |\/b 



70 



=P4>Ji 

(27/3 (3-1 + 2 7 2 ) - # 2 ) J 2 



16|p 2 |VTFT V Q l+7 2 7~ u ' 16|p 2 | v /b T |(l + 7 2 ) V \ ol 
This concludes the classical symmetry reduction of canonical gravity non-minimally coupled to fermions. 



(21) 



E. Parity behavior 



Because we are mainly concerned about the role of parity in loop quantum cosmology, we end this section on the 
classical equations with a discussion of parity invariance. As pointed out in [l6j . parity invariance in loop quantum 
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gravity is not guaranteed. The Ashtekar connection is a sum of a parity-even and a parity-odd term and thus does not 
have a straightforward parity behavior. This already occurs in the absence of fcrmions and torsion, but is aggravated 
by the parity-mixing terms of torsion contributions due to a fermion current (see (|A6|) , noting that J\ is even and Jo 
is odd). Classically, one can explicitly split these contributions, which essentially amounts to replacing the Ashtekar 
connection with extrinsic curvature. However, a complete splitting requires equations of motion to be used, which will 
not be possible after quantization. It is then not guaranteed that quantum corrections due to the loop quantization 
will preserve parity even in vacuum or in the absence of parity-violating matter. 

The model introduced here provides a clear view on parity in the classical theory as well as after quantization, as 
we will see below. One key property is that the canonical variables (jlll) we are led to do, in hindsight, perform the 
splitting into torsion-free and torsion components without using equations of motion. Thus, in the new variables every 
single term in the Hamiltonian constraint (|21[) has a clear and simple behavior under parity: Among the gravitational 
variables, only c\ and p 1 change sign under parity (reversing orientation) while the rest remains unchanged. (Since 
changing the sign of p 2 in an LRS model implies a reflection of both directions related by the rotational symmetry, 
it is equivalent to a triad rotation and thus mere gauge.) This is accompanied by the usual parity transformation 
of the fermions present, which implies that J7i is parity invariant while j7o changes sign as these are space and time 
components of an axial vector. In particular, it is immediately clear from (|2ip that the Hamiltonian constraint is 
parity invariant for free fermions. Parity violation will only result if suitable interactions are introduced to the model, 
which can easily be done by adding e.g. y/— dct gV^J^ with the vector current V 1 to the action. We will avail 
ourselves of this possibility in what follows to understand the role of parity in the loop quantized model. 

III. QUANTIZATION OF THE BIANCHI I LRS MODEL 

Loop quantum cosmology allows one to complete many of the constructions of full loop quantum gravity in simplified 
and explicit forms, which then provides indications toward the physical implications of the theory. In this section, 
we provide a self-contained description of anisotropic models with an emphasis on the effects of fcrmions, torsion and 
parity. 

A. Quantum Kinematics 

We start with basic variables according to the Poisson structure (|T7|) . As in any loop quantization, states in the 
connection representation are constructed by taking exponentials 

exp(/iiCiAj n) £ SU(2) , exp(/i 2 a 2 A 4 2 T J ) £ SU(2) , exp(zfc0) £ U(l) for all m £ K, k £ Z, A} £ 50(3) 

as they arise in holonomies. Using holonomies in the general setting is important for a background independent 
basic algebra of variables. This crucial feature is then reflected also in symmetric models based on exponentials of 
connection components. The parameters pi can take any real value, corresponding to evaluating holonomies along 
straight edges (tangential to Xf) of arbitrary length. The variable (f>, on the other hand, was introduced as a periodic 
angle in ©, such that only strictly periodic functions exp(ik(f>) with k £ Z are allowed. This unphysical degree of 
freedom, which we were led to introduce due to the presence of torsion, will be removed after solving the Gauss 
constraint. 

Matrix elements of the exponentials in (|22|) form a C*-algebra of (almost) periodic functions, as seen from ([5]). Any 
function generated by this set can be written as 

g(c 1 ,C2,<f>) = ^ £/*i,P2,ife ex P (l^ici + i«>2C2 +ik<(>) , (23) 

with coefficients £ C, where the sum is over finitely many fj>x, /j>2 £ R and k £ Z. Note that while g(ci, £2, </>) 

is almost periodic in c\ and c 2 , it is exactly periodic in 4>. This provides a complete set of continuous functions 
on Reohr x RBohr x S 1 , where RBohr is the Bohr compactification of the real line. (By definition, RBohr is the 
compactification of R such that the set of all continuous functions on it is just the set of almost periodic functions. 
See e.g. [32[ for a recent discussion of further properties.) All spaces in the product are compact Abelian groups and 
carry a unique normalized Haar measure d/i(c) in the case of R.Bohr ; where 

//(ewe) '.= **±f m*c, 
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and d(f) for S . 

By Cauchy completion, we obtain the Hilbert space as a tensor product H grav = ^fohr ® ^s 1 with the Hilbert 
spaces HBohr = L 2 (RBohr, d/i(c)) and Wgi = L 2 (S rl , d^) of square integrable functions on the Bohr compactification of 
the real line and the circle, respectively. Orthonormal bases for these spaces are given by (c\fi) — exp(z/xc/2), fj, € K, 
and (0| fc) = exp(ifc^), fc G Z, respectively, with 

0V> = <w . = ^,fe' • ( 24 ) 

The configuration variables act in the obvious manner: For all g\ and <? 2 of the form (|23|) . we have 

(?i52) (ci,C2,4>) = 9i(ci,c 2 ,4>)g2(ci,C2,(f>) (25) 
and the momentum operators are represented by 

, p 2 = -4^ - (26) 

where £p = kH. (The densitized triad in general is quantized via fluxes, i.e. 2-dimensional integrations over surfaces. 
In a homogeneous context, however, this is not required and densitized triad components can directly be promoted to 
operators. This simple representation exists only due to our use of variables; had we used c 2 instead of c 2 , the operator 
p 2 and thus the volume operator would have been more complicated.) Common eigenstates of all triad operators p 1 
are 



\Hi,H2,k) := |a*i> ® |/i 2 ) <S> \k) , (27) 

with 

p 1 \lH,H2,k) = -^Y^- lAH,M2,fe) , P 2 W,^2,k) = ^-^- \fii,H2,k) and p</,|/ii, At 2 , fc) = 7£pA;|/ii, /i 2 , fc) . (28) 
From triad operators we construct the volume operator: 

V\fii,fi2,k) = IflVWl \vi,fJ-2,k ) = — ^|/Li 2 |vM |Mi,M2,fc) • (29) 



4V2 

The full Hilbert space is a further tensor product of H gra v with the fermionic Hilbert space H fermion- We represent the 
latter as the space of functions f(Q a ) of four independent half-densitized Grassmann- valued variables 6 a , a = 1, ... ,4, 
for the four components contained in the fermion fields £ and \ m this order. The fermionic momenta ir^ = and 
ir x = — then give rise to components Q a which are represented as hd/dO a . In particular, for the axial current 
components J° = £t£ — x^X an d J7i = £ t(J i£ + X fiX we have operators 

oBi oB 2 063 064 

oB 2 aBi 0B4 aB 3 

(The component j7o is subject to ordering ambiguities which we can ignore here.) 

The currents are easy to diagonalize: Each 2-spinor copy has eigenstates of g§^©i + g§7@2 given by /o(B) = 1 
and /°(0) = BiB 2 of eigenvalue zero, /±(B) = Bi ± B 2 of eigenvalue ±1. The tensor product of both 2-spinor 
copies £ and \ then gives eigenstates of eigenvalues zero, ±fi and ±2S for J\. The time component j° has the same 
eigenstates. 

A general state in TL = 7i gra v ® ^fermion can then be written in a form using fermion dependent coefficient functions 
in the triad eigenbasis (|27|) : 

I s ) = ^2 s w , M2 ,fc(G)|/xi,/i 2 ,A;) . (32) 

One can define the coefficients Sui,u 2 ,fe(®) f° r au " values of /ii,/i 2 G K and fc € Z in this way. However, gauge invariance 
implies that the state must be invariant under changing the sign of /u 2 because this corresponds to a triad rotation 
(without changing orientation). Thus, we require s All , A12 .fe(B) = s Ml ,-^ 2 .fc(B). 
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The remaining sign freedom, sgn/xi, is physical and crucial because it determines the relative orientation of the 
triad. Thus, we have a simple action 

S Ml ,^,k(&l, ©2, ©3, ©4) — > s_ AIl!A12! fc(6 3 , 84, 61,82) (33) 

of the parity operator II on states. For the fermion dependence, we have represented the parity action ITO = 7°^ for 
Dirac spinors by switching the fermion values 8 Q corresponding to £ and Xi respectively. This implies 

njoft- 1 - -Jo , njxrr- 1 = j x . (34) 

For gravitational operators, a direct calculation shows 

ft^fr 1 = -p 1 , tiffir 1 =f (35) 

flexp(i/uici/2)n _1 = exp(— i/iiCi/2) , ITexp(i/i2C2/2)n _1 = cxp(i^ 2 C2/2) (36) 

as required. 

Finally, we can directly solve the Gauss constraint which requires p,p = h'yuJi and thus allows us to eliminate k as 
an independent quantity. Using the spectra of the operators already determined, this provides solutions with either 
k = or A; = ±1. In the second case, there is a non- vanishing value of the spatial axial current J\ of size ±2ft. The 
values ±h for the fermion current, which do exist as eigenvalues, are ruled out because they do not correspond to 
integer k. Both 2-spinors present must thus have the same or opposite J\ -eigenvalues, which allows them to be parity 
eigenstates. The parity behavior of the full state according to (|33p . however, is determined by the /xi-dependence, 
which required the dynamics of quantum gravity coupling the triad to fermions. 

The allowed values for the current are only microscopic and may not seem of interest to describe a macroscopic 
universe of large matter content; they all vanish in the classical limit h —> 0. Nevertheless, this provides an interesting 
model where one can study the effects of fermions and parity in loop quantum gravity. Physically, it is also clear why 
the matter contribution can only be microscopic: As always in homogeneous quantum cosmological models, each held 
component is reduced to a single degree of freedom for all of space. For the fermion, this allows only one excitation 
per component due to Pauli's principle. Unlike with scalar matter, one cannot simply make the matter content large 
by choosing a high "occupation" such as a large momentum of the scalar. Significant fermionic matter can only 
be included by adding more independent spinor fields, or by introducing inhomogeneity which provides independent 
field values at different points (represented by fermions at different vertices of a spin network state in loop quantum 
gravity). Rather than being a limitation, we consider this as an important physical property of quantum cosmology 
in the presence of realistic fermionic matter. 

B. Quantum Dynamics: The Hamiltonian Constraint 

A useful feature of the torsion-free Bianchi I model is that the Lorentzian Hamiltonian constraint is related to the 
Euclidean part simply by H = —j~ 2 H( E ' ) thanks to K^Kh oc F^ b tijk, making use of homogeneity as well as the fact 
that the spin connection vanishes. This has been used in almost all investigations of loop quantum cosmology so far. 
If this relation is not used, one can still quantize the Lorentzian constraint following techniques of the full theory 33]. 
This results in a more complicated constraint operator but without crucial differences. 

However, in the presence of torsion, such a simple relationship can be obtained only after splitting the torsion 
contribution from the spin connection as shown in (|A9[) . which is now to be quantized: even for the Bianchi I model, 
r* is no longer zero due to torsion. Fortunately, torsion contributions to r*, namely C l a in (|A6[) . are completely 
determined by second class constraints. They can thus be split off and quantized separately together with the matter 
terms. For the Bianchi I LRS model, one can use a further key simplification which, as pointed out above, allows us 
to project out torsion contributions without directly computing them. All we need to do is use the new variable C2 
instead of ci. The resulting contribution to the gravitational Hamiltonian constraint is the same as the torsion-free 
one and thus can be quantized in the same way. 

Mimicking the steps done in the full theory [TBI. [33|. one writes curvature components T % ab as a product of (point) 
holonomies hi = cos(i<5fCf) + 2A}TjSin(i5/C/) forming a closed loop, whose "edge lengths" are denoted as <5i and <5 2 
for the two independent directions. Moreover, using 
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relevant products of triad components, including their inverse powers, are reduced to a Poisson bracket of the general 
form hilhj 1 , V} where V is the spatial volume and hi again a holonomy. This allows one to write an operator in 
compact form, which corresponds to a densely defined operator in the full theory: 

" -iSfKf £ ^ (^A-^'M^]) ■ (38) 

' UK 

We can now compute the product of holonomies and take the trace explicitly, using the basic properties of Pauli 
matrices. We do this directly for LRS variables with only two independent holonomies such that 82 =83. Moreover, 
the sign factor is now solely determined by sgnp 1 since p 2 p 3 cannot be negative. This results in [2(| 

H g = - j3 K £2 § $2 (^2sin(i(5iCi)cos(i(5ic 1 )sin(i(5 2 c2)cos(i(52C2) (sin(|5 2 c 2 )y cos(i£ 2 c 2 ) - cos(±5 2 c 2 )V sin(i<5 2 c 2 

+ sm 2 (i<5 2 c2) cos 2 (i<5 2 c 2 ) f sin(i£iCi)Fcos(|#iCi) — cos(|<5iCi)F sin(^8iCi)jj . (39) 

Because we have implicitly eliminated the torsion contributions from holonomies by our choice of basic variables, we 
can directly use this expression as it is known from torsion-free models. The torsion contribution will then be added 
to the constraint operator via the fermion current. 

We emphasize that the meaning and form of the parameters 8\ and 82 cannot be fully elucidated purely in homo- 
geneous models. In the absence so far of a derivation from a full, inhomogeneous constraint (which itself is currently 
subject to changes in its general form depending on ongoing developments) it appears best to refrain from specific, 
heuristic arguments as to what values they may take. (For instance, there is currently no firm basis for a relation of 
those parameters to an eigenvalue of the area operator of the full theory, as initially proposed in [3J].) We therefore 
follow a more general route which allows whole classes of these parameters, and confine attention to effects which 
are insensitive to the specific form. To us, this seems most advisable given that it is not just the numerical values of 
these parameters but even their possible functional dependence on basic variables which remains open; see Sec. IIII CI 
for further discussions. 

In order to quantize the matter Hamiltonian, we must in particular quantize the inverse volume ^/p 2 \/\p 1 \- Here, 
we use the standard procedure [l5[ , first writing 

1 _ sgndct(e^) ■ a j k_ 36 scrnt ictrcMc abc c i\A l V 1/3 \\A j V 1/3 \iA k V 1/3 \ 
V |dct(£?)| ~ 6|det(£f) ' jtWc "7 3 K 3 g l J ^k\A a ,V i\A b , V i\A c , V j 



based on (|37|) . which is then quantized to 

l\ 144isgn(p 1 p 2 p 3 ) ^ IJK 



^ C " " tr [hi [hj 1 , V 1/3 \ hj [hj\ V^ 3 ]h K [h^ 1 , V 1 ' 3 ' 

UK 



sin(i(5iCi)T> 1/3 cos(i<5iCi) - cos(±6iCi)T> 1/3 sin(±(5iCi) 



32 • SlsgnQj 1 ) 
~1 3 1%8 X 8 2 

(sin(|5 2 c 2 )F 1/3 cos(i<5 2 5 2 ) - cos(±8 2 c 2 )V 1/3 sin(±<5 2 5 2 )) . (40) 



The action of this operator as well as the Hamiltonian constraint is easily computed using the action of sin(i<5iCi) 
and cos(MiCi) on the triad eigenstates, 

cos(i6iCi)|/ii,^2,fc) = ^(|mi +£i,M2,fc) + iMi - Si,H2,k)) 

sin(i5ici)|//i,^ 2 ,fc) = --i(\ni+6i,H2,h) -\fj,i-6i,(jQ,k)) , (41) 

and the volume operator (j29|) . From matrix elements of the Hamiltonian constraint one can then write the constraint 
equation (Hq + -ffmatter)|s) = as a difference equation for coefficients s Ml , M2 ,fc(0) of the state in the triad represen- 
tation. We do this immediately on states solving the Gauss constraint which determines k in terms of the action of 
J\. Dropping the label k on those states, we have 

2(|M2 + 3S 2 \ - \fi 2 + 8 2 \) +2(5ir /2 s Ml+2 5 1 , At2 +25 2 (e)- \m -2<5 1 | 1 / 2 s Ml _2 5l , M2+ 2i 2 (e) 
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+2(|/i 2 - 5 2 \ - |M2 ~ 3<5 2 |) (| M i - 25 l \ 1 ' 2 s^2S 1 ,^-2S 2 (Q) - |mi + 2<5i j 1/2 s Ml+2<5l _ 2(52 (O)) 
+Oi + 1 1/2 - Imi - <5i| 1/2 ) (lM2 + 4£ 2 |s Ml , M2+453 (e) - 2|^ 2 |s Ml , M2 (6) + |^2 - 4<5 2 |s W)M2 _45 2 (e)) 



= ^M 1/3 M 1/3 (lm + ^i| 1/6 - iMi - ^i| 1/6 )(lM2 + <5 2 | 1/3 - |M2 - ^| 1/3 ) 2 
16 

27/? / g 7 2 \ A 2 \ j? 
1 + 7 2 V a ) 1 +7 2 / h 2 ' " '" V" ' 1 + 7 2 



X |(l+V-^3- ^-^ ^^ +T ^ § , tt ,(e). (42) 



This equation is based on a non-symmetric constraint operator because in (|39p we ordered all holonomy factors 
to the left and kept the commutator terms with the volume operator to the right. It is sometimes useful to have a 
symmetric ordering, although this is not strictly required for constraints. (But it is required by some methods to 
derive the physical Hilbert space.) There is only one way to order the constraint symmetrically, namely by introducing 
\{H + iF). Other possibilities have been suggested, such as splitting the sines and cosines and writing some to the 
left, others to the right of the commutator term. They are, for instance, useful to prove self-adjointness [351 ] . However, 
this corresponds to splitting the holonomy product hihjhj 1 ^ 1 into different factors, which cannot be done in a 
general setting where there would rather be a single holonomy h a around a closed loop a. The direct symmetrization, 
on the other hand, is always possible and in our case results in a difference equation 

2 ((|/i 2 +3S 2 \- | M2 + hDlm + 2<5 1 | 1 / 2 + (| M2 + S 2 \ - \fx 2 - 5 2 |)|mi| 1/2 ) s^ +25l ,^ +252 (e) 
-2 ((|/Z2 + 36 2 \ - | M2 + S 2 \)\ fH - 25 1 \ 1 ' 2 + (| M2 + 6 2 \ - |/x 2 - £ 2 |)|mi| 1/2 ) s w - 25l>(U2+2& (e) 
+2 ((| M2 - 5 2 \ - |M2 - 3<5 2 |)|Mi ~ 2<5i| 1/2 + (l/fc + S 2 \ - | M2 - ^ 2 |)Imi| 1/2 ) V-a*,/*-a«,(e) 
-2 ((|M2 - S 2 \ - | M2 ~ 35 2 |)|/ii + 2,5 1 | 1 /2 + ( | M2 + 5a j _ | M2 _ 5 2 |)| Ml |i/2^ fl|tl+Mlf ^_ Ma (0)) 

+ (iMi + ^i| 1/2 - |A*x ~ <5i| 1/2 ) ((H + |/i 2 + 4(Si|)a Ml , M+ 4« 2 (O) - 4|/x 2 |s w , M2 (e) + (| M2 | + |/* 2 - 4*i|)a w ,»-4*i(Q)) 

= f ImiI 1/3 Im 2 | 1/3 (Imi + *i| 1/6 - |M! - ^il 1/6 )(lM2 + s 2 \^ - | M2 - 6 2 \^f 

2 7/ 3 / 7 , 2 \ 2 \ J? , q fl ^2 , 7 \ 



Lattice refinement 



So far, we have left the increments <5i and 6 2 unspecified. It is clear that as constants they would not influence the 
recurrence behavior of the difference equation, although specific solutions certainly depend on their values. However, 
in general and S 2 may not be constant but be functions of [i\ and fi 2 ; this captures the way in which the discrete 
structure of a state underlying spatial expansion and contraction in loop quantum gravity is being refined dynamically 
[26L l36l [37l| : at larger /ij, an increment of the total size by a Planck-scale amount has a weaker relative influence on 
the geometry. As a consequence, 5i decrease with increasing spatial extensions. This can also be seen from more 
direct considerations of holonomies in inhomogeneous states and how they appear in Hamiltonian constraint operators. 
Since this involves the dynamical relation between models and a full non-symmetric theory, the precise behavior of 
lattice refinement has not been completely determined. However, consequences of different behaviors can be explored 
in several models. Sometimes, this is already quite restrictive even though it is impossible to derive a unique form of 
lattice refinement based solely on homogeneous models. 

Non-trivial functions, such as power laws, have a much stronger influence than constants because they make the 
difference equation non-equidistant. Solutions are then more difficult to analyze and find, even numerically (but see 
[38l [39(1 ^ . Only in the special cases where 5\ cx /i^ 1 and S 2 oc /i^ 2 can the equation be mapped to an equidistant 
one by a redefinition of independent variables. However, such cases have been ruled out [401 ] because they do not 
provide the correct semiclassical behavior near a horizon of Schwarzschild black holes, whose interior is treated as 
a homogeneous Kantowski-Sachs model. (The analysis in [40( uses corrections to classical equations due to the use 
of holonomies in the loop quantization, but ignores other effects such as quantum back-reaction [4ll |42|. This type 
of phenomenological equations may not capture correctly the behavior of strong quantum regimes such as the black 
hole singularity. However, if these equations do not provide the correct semiclassical behavior in classical regimes, 
this cannot be corrected by the inclusion of quantum back-reaction. The fact that some refinement schemes are ruled 
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out is thus a reliable feature.) In general, one has to expect functions of the form <5i(/xi, and (^(mI) M2) with a 
non-trivial dependence on both arguments (which may not be of power-law form). 

As we will see below, a discussion of fundamental singularity resolution only involves the recurrence near /ii = 0. 
This is, fortunately, insensitive to the particular refinement scheme and thus presents a result of much wider generality 
than anything which applies at larger volume where the specific refinement can be crucial. 

IV. COSMOLOGICAL IMPLICATIONS 

It follows immediately from the difference equation (|42f or (l43|) that it is parity invariant since all its terms change 
sign under (|33|) . Thus, if s Ml!M2 (81, 82, 63, 84) is a solution, so is s- lll ^ 2 (83, 84, 81, 82). In particular, any solution 
can be written as a combination of even and odd solutions 5^^(81,82,83,84) ± S-^^ (83, 84, 8i, 82). This is 
no longer the case if we had matter interactions violating parity, such as a term proportional to VoJo- In this case, no 
parity-even or odd solutions would exist. Wave functions for fix > generically differ from their form for /ii < 0, even 
though those values are deterministically related via the difference equation. At this stage, the precise form of parity 
violations in the matter system is crucial to determine the behavior of the wave function near the classical singularity 
at fii = 0. 

To complete the construction, one would solve the difference equation and determine a physical inner product on 
the solution space. Ideally, one could then compute the behavior of observables of the system and derive detailed 
cosmological scenarios including the role of quantum effects. Unfortunately, such complete descriptions at an exact 
level are possible only in rare, specific models. While such models are instructive mathematically, conclusions drawn 
are difficult to interpret because one could not be certain about the robustness of results: If specific results are 
available only in a few special models where exact mathematical solutions in the physical Hilbert space can be found, 
there is no guarantee that they are not just the very result only of demanding this high mathematical control. 

In this context, an aspect of particular interest is the fact that most models of loop quantum cosmology where 
physical Hilbert spaces have been constructed explicitly [l3|, |43|, |44[ specifically assume parity invariance in some form 
and make use of the corresponding restriction of states when parity is considered as a large gauge transformation. As 
we have seen here, physical states of quantum cosmology are neither even nor odd in triad reflections if parity violating 
matter is present. It may thus be misleading to treat parity as a large gauge transformation even in cases where matter 
preserves parity. Results based on this assumption may be spurious, and one has to re-analyze the constructions of 
physical Hilbert spaces without the assumption of parity invariant states. Fortunately, the intuitive pictures of 
bounces which have sometimes been derived from physical observables are insensitive to the specific construction of 
the physical Hilber t sp ace: They can be derived analytically in a representation independent formalism based on 
effective equations jlOl . |45| . Then, the assumption of parity as a large gauge transformation is not necessary, and it 
can be dropped without affecting the bounce result. 

At a fundamental level, singularity resolution is also insensitive to the physical Hilbert space construction and 
can directly be determined using the difference equation PS|) or (pt3")) . (Here, it is important that all solutions are 
non-singular, which then also includes physical ones.) In general, coefficients of a difference equation of the type 
obtained in loop quantum cosmology may vanish and prevent certain values of s Ml /J2 from being determined in a 
recurrence starting from initial values. This happens for the non-symmetric equation (I42|) where none of the values 
'00,^2 — right at the classical singularity — is determined by initial values because their coefficients in the difference 
equation vanish. (The corresponding states |0, ^2) are mantic @.) However, for the difference equations realized such 
undetermined values, if they arise, drop out completely of the recurrence. In particular, even though values for /xi = 
remain undetermined by initial values in the non-symmetrized version of the equation, coefficients at jix < follow 
deterministically from coefficients at fii > 0. 

In parity preserving models the wave function s Ml . M2 for fii < could simply be the mirror image of its cousin at 
fi± > 0, and it had to be symmetric if parity is considered a large gauge transformation. However, if there is parity 
violation, the transition through fii = constitutes true evolution since values at fXi < must now differ from the 
mirror image at /ii > 0. The wave function at \i\ < cannot be determined simply by reflection, but it has to be 
derived by local evolution through all intermediate values of In this case, the region of fii < can by no means be 
removed from considerations but must be considered as a physical domain on equal footing with that at [i\ > 0. In 
particular, the orientation-reversing big bang transition thus becomes physical and cannot be argued away as a large 
gauge transformation. 

For both forms of difference equations derived here, there are consistency conditions arising due to vanishing 
coefficients around fi\ = 0, analogous to dynamical initial conditions [111 Il2j . If we evaluate any of the difference 
equations at \i\ = 0, matter terms drop out and we obtain the universal relation 

(|/i 2 + 3<5 2 | - |/-*2 + S 2 \) S2S 1 ,fj, 2 +2S 2 - (1^2 - S?\ - |M2 - 3(5 2 |) S2S 1 , f i 2 -2S 2 
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= (|/i2 + 35 2 | - |M2 + S 2 \) S_25 1)/i2 +2« a ~ (|M2 ~ fe| - 1^2 - 3<5 2 |) S_ 2 5 1 , AI2 -2<5 2 (44) 

valid for all /i 2 . In particular, at /i 2 = 2<5 2 we have s 2< 5 1 ,4<5 2 = s^2S 1 as 2 - At odd integer multiples of /i 2 = 2<5 2 , we 
obtain a recurrence relation which requires s 2( 5 lj2 ( 2n+ i),5 2 = s_25 1 ,2(2n+i)(5 2 f° r a ^ integer n. 

There are thus reflection symmetry conditions which directly follow from the dynamical law even in the presence 
of parity- violating terms. (This symmetry has been observed first in the vacuum case 46].) However, evolution away 
from fix = ±1 depends on whether \x\ is positive or negative if parity is not preserved. Thus, the wave function is not 
mirror symmetric even though the dynamical initial condition closely ties the values s±2S u ti 2 to each other. 



V. CONCLUSIONS 



We have introduced fermions into the framework of loop quantum cosmology which gave rise to several non-trivial 
changes due to the presence of torsion and potential parity non-invariance. We have observed several key features 
which have a bearing on cosmological scenarios and which do not arise for bosonic matter such as scalar fields as 
they are used commonly in cosmological models. First, the amount of matter is limited for each fermionic degree of 
freedom due to the exclusion principle. Thus, large matter contents as they are sometimes used to bring a quantum 
cosmological model into a semiclassical regime where it may bounce more easily cannot straightforwardly be achieved. 
The only possibilities are to allow many copies of independent fermions or inhomogeneity where fermionic components 
at different points will be independent. Physically, both possibilities are quite different from having a single bosonic 
field of high occupation. The methods used here may also be of interest for a supersymmetric version of loop quantum 
cosmology along, e.g., the lines of [13] (see also [HJ). Fermions in quantum cosmology also play a role for decoherence 

This shows that it is crucial to consider the small- volume regime of a quantum cosmological model which cannot be 
avoided in the absence of much matter energy. Here, the recurrence scheme of an underlying difference equation of loop 
quantum cosmology becomes essential to determine whether the model is singular or not. As we showed, the singularity 
resolution mechanism of loop quantum cosmology [6j remains unchanged under the inclusion of fermionic matter even 
if it violates parity. At the same time, the model we used allows us to show that in its realm parity violations can 
only arise due to matter interactions, not due to pure gravity. In other models or the full theory, this situation may 
be different because the basic objects quantized, in particular holonomies, do not transform straightforwardly under 
parity. The model introduced here thus also serves the purpose of providing one example where parity invariance of 
pure gravity can be demonstrated after a loop quantization. 

If one introduces parity-violating interactions, wave functions cannot be mirror symmetric. Then, the branches 
at the two opposite orientations of triads are independent of each other, and joined through degenerate geometries 
by the dynamics of loop quantum cosmology. The big bang transition now becomes a non-trivial event where space 
turned its inside out in a quantum process which in general cannot be described by an intuitive geometrical picture 
such as a simple bounce. 
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APPENDIX A: THE FULL CONSTRAINTS 



To set the notations, the basic configuration variables in a Lagrangian formulation of fermionic field theory are the 
Dirac bi-spinor = (-0, rf) T and its complex conjugate in <3> = (\E'*) T 7 with j a being the Minkowski signature Dirac 
matrices. We note that tp and 77 transform with density weight zero and are spinors according to the fundamental 
representation of SL (2, C). Then the non-minimal coupling of gravity to fermions can be expressed by the total action 
composed of the gravitational contribution Sg and the matter contribution Sf resulting from the fermion field: 

S[e,w,V] = S G [e,oj}+ S F [e,oj,t>} (Al) 
1 f d 4 x \e\e^e v J P IJ KL F,J CL ( U ) + \i [ d 4 z |e| (% 7 'e? ( 1 - -75) - V^" ( 1 - - 75 ) /e^ 



lD7rG J m * 2 J M \ \ a J \ a 
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Here I.J,... — 0, 1, 2, 3 denote the internal Lorentz indices and fi, v, . . . — 0, 1, 2, 3 space-time indices, and a € M is 
the parameter for non-minimal coupling. 

We have expressed the gravitational action in terms of the tetrad field ej , where e is its determinant and the 
inverse, using the Hoist action [50]. It presents a first order formulation where the Lorentz connection, denoted by 

uj 1 ^ 1 and with curvature F^ v l {uj) = 2d[ fJi ui I i/ ^ + , is treated as a variable independent of the tetrad before 

equations of motion are imposed. In the Hoist action, we have 

P IJ kl = S [ X ] -Z^ (A2) 



with inverse 

-.2 / ... _. 1 

t IJ 



IJ — „,2 , i °J °J 



where 7 is the Barbero-Immirzi parameter [EH |52| . The connection appears also in the matter part via the covariant 
derivative V M of Dirac spinors defined by 

V M . = 9 M + ^ J 7[/7j] , [V M ,V„] = ^ 7 7[/7j] (A3) 

in terms of Dirac matrices 7/ (which will always carry an index such that no confusion with the Barbero-Immirzi 
parameter should arise). 

For a — > 00 we have minimal coupling from the viewpoint of the Hoist action, while (1 = 7 corresponds to minimal 
coupling from the viewpoint of Einstein-Cartan theory [23|. In this article, we allow all possible real values of a as 
introduced in [23| (see also [Hj]). We emphasize that we have parity invariance for all real a, even though some torsion 
components such as those written below consist of contributions of different parity behavior. 

A canonical analysis of the action yields first and second class constraints. To summarize the result of this analysis 
[l6| . we use k — 8irG, tj — ~^<Jj in terms of Pauli matrices, the axial fermion current components 

J 1 = V>VV + T7V77 , J° = - rfv (A4) 

and the parameters 

= 1-1 and /3 = 7 + - (A5) 
a a 

which are useful as a shortcut. For a consistent loop quantization, the half-densitized £ := tfqip instead of (and 
X := tfqr) instead of 77) is required to be the classical canonical variable for fermions and tt^ = — i$ is the 
conjugate momentum for £, using the spatial metric q a b and its determinant q. 
Upon solving the second class constraints, which provides the expression 

°i=4(i?^)(^«^-^) (A6) 

for the torsion contribution to the spin connection V a , the usual first class constraints remain: the Gauss constraint 

G[A] := / d 3 xA* (V b P* - \y/qJi) = f d 3 a;A l (V b P t b - tt^ - 7t x t iX ) , (A7) 

the diffeomorphism constraint 

D[N a ] := £d 3 a; iV^Pj'S^-^^ (A8) 
and the Hamiltonian constraint (modulo Gauss constraint) 

HtoUN] = jf cl 3 x N (^|i?P/ (e« fe ^ 6 - 2( 7 2 + 1)4,**) + 1 ^^aJ l + ^Pfy 
U~*» c „S» „ _ . , ^ l!*!^ p. J. *o„ pb _ 37^ / 2 ^ 7^ 



(27^(3- 2 + ^)-*)^) (A9) 



16 



expressed in terms of the densitized axial fermion current components 

^=eVe+xV x , j° = ?t-x*x- (aio) 

The Hamiltonian constraint has been written in a useful form after splitting the spin connection T* = T* + C l a with 
the torsion-free connection T % a . As usually, N and N a are, respectively, the lapse function and shift vector used to 
foliate the space-time manifold M. 

As basic gravitational variables we use a canonical pair given by the densitized triad ff = Ef/'yn together with 
the connection 

7K 



A ■■= A i+j^< j0 =K+ci+ 7 Ki (ah) 



ei = & a + 5:4.7° = ( \ Jk ^j _ eU o ) (A12) 



where 

C l — C 1 -L _ 

4a" a " 4(1 + 7 2 ) V7 

whose curvature and covariant derivative we denote as T^ b and T>, respectively. The J°-term in the connection A l a , 
compared to the Ashtekar-Barbero connection, is required for a formulation in terms of half-densitized fermions. 
Solutions of the second class constraints tell us only what the torsion contribution C l a to the spin connection is. To 
know the torsion contribution to extrinsic curvature K l al and thus to the Ashtekar-Barbero connection, one has to 
partially solve equations of motion. Doing so [IB], the J°-term is cancelled and we can write 

A\=fl +1 Kl + ^ kl e k J l . (A13) 
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